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We numerically study energy spectra and localization properties of the double exchange model 

^— «, , at irrational filling factor. To obtain variational ground state, we use a mumerical technique in mo- 

("—^ ' mentum space by "embedded" boundary condition which has no finite size effect a priori. Although 

r\] \ the Hamiltonian has translation invariance, the ground state spontaneously exhibits a self-similarity. 

Scaling and multi-fractal analysis for the wave functions are performed and the scaling indices a's 

(^ ' are obtained. The energy spectrum is found to be a singular continuous, so-called the Cantor set 

with zero Lebesque measure. 
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The study of the double-exchange (DE) model ||l| has a long history. The purposes are to understand the elec- 
tronic P] and transport properties H of a class of the transition metal oxides, mainly the perovskite nianganites, 
which are known to show the colossal magnetoresistance (CMR) 0]. Recent experimental and theoretical studies 
^ ■ have revealed that it is a strongly-correlated quantum liquid in which many degrees of freedom are participating, 
I ] such as orbitals, lattice distorsion, electron correlation, dimensionality, the quantum nature of the spin, etc pM]. The 
^ . current researches are in progress towards the investigation of the systems with many degrees of freedom. A variety 

a of "double-exchange model" have been introduced and investigated from many different viewpoints and approaches. 

. . On the other hand, the minimum model of the DE mechanism has natures still to be clarified even before introduction 

of the many degrees of freedom. (The model is composed of the itinerant Cg electron and the localized t2g spin coupled 
with the Hund coupling each other.) One is phase separation p|-p^. They has been found not only in the minimum 
S-H ■ model also in those with direct exchange couplings between t2g spins |9| with orbital degree of freedom, and/or with 
i-rt i electron correlation W^. Another is the spin-induced Peierls instability |l^JT3]. The ground state is characterized 
r^ ■ by a charge gap and a modulated spin structure. The charge gap is induced by a commensurate modulation of the 
O ' transfer integrals. The spin state is distorted with a period commensurate with the Fermi momentum to open the 
, ^, , charge gap, thus stabilizing the system. This is isomorphic to the Peierls instability ||l^ which has been discussed 
for a long time. The only difference is that the spin degree of freedom is distorted, instead of the lattice degree of 
ff^ , freedom. 

^ I In this letter, we numerically study the electronic structure of a one dimensional DE model at an irrational filling 

f^ I factor. Although the Hamiltonian has a translation invariance, the energy spectrum spontaneously exhibits a self 
^^ ■ similarity and forms the Cantor set. The scaling analysis for the band width and the wave functions are performed. 
2^ [ We estimate the exponents a's which characterize the localization properties [|T^. To study the fine structures, usual 
f—N . boundary conditions, such as periodic, twisted or open one, are useless. We use an alternative numerical technique in 
f~*) ■ the momentum space with "embedded" boundary condition. This has no finite size effect a priori. The localization 
f^ , properties and the similarity with the quasi-periodic system are discussed. The electronic and transport properties 
^^ ■ are found to be completely different from those of the DE model in the literatures. 
Cu ' The DE model we employ is 

g: 

Y . Li/2 

• • , where Ci^a- is a fermion (cg electron) annihilation operator at site i, Si are localized (t2g) spins which are treated as 

K*" • classical vectors directed along {9i,(j)i) in the spherical coordinates. Moreover, J>0 is the direct exchange couphng 

k>( ' strength between t2g spins. We fix the parameter to be JS'^/t=l throughout this paper. L (even) is the total number 

j_j of the sites. The transfer integral |l^ is given by 
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The pure one-dimensional version of the model M), i.e. in absence of the second term, all the sub-bands are disper- 
sionless |L4| . We introduced the second term only as a perturbation to remove the singular behavior. The singular 
behavior can be removed by any perturbation and the form of the second term is not special to produce the band 
width. 

We study incommensurate structure as a systematic limit of commensurate ones. The difficulties we meet are the 
followings: (i) The spin structure of the ground state is expected to be incommensurate, which are very sensitive to 



boundary conditions. When we impose periodic, open, or twisted one, the incommensurate structure is restricted to 
the corresponding space. The true grournd state for a finite size does not necessarily exist in the space, (ii) The 
periodic boundary condition forces to aUgn the spins periodically and spoils the fine structure. One might expect 
that the incommensurate structure can be asymptotically embedded into the commensurate ones in a sufRciently 
large system. However, the finite size correction is roughly estimated for 1/i, which is not small for our purpose. 
The free boundary condition often induces a boundary state and provides the similar problem, (iii) We need a large 
denominator of the filling for the scaling analysis, (iv) The ground state has a continuous infinite degeneracy. To 
avoid these difficulties, we use "embedded" boundary condition to the spin part. The boudnary condition is similar to 
Ref. jlj]. We set the oriented relative angle between Sl and Sl+i is the same as that between iS'o=(0, 0) and ^i. (The 
spins of i>L+l are successively redefined from those of l<i<L+l with maintaining the oriented relative angles.) We 
impose the periodic boundary condition for the eg electron, i.e. cl+i=ci. This enable us to treat the non-translation 
invariant spin states under the Fourier transformation for the Cg electron. No finite size correction due to (i) and (ii) 
exists a priori. We further assume the spin-induced Peierls conjecture |l3|JI^ ]. Finally, at filling x=p/{2q) the kinetic 
term of the Hamiltonian in the momentum space reads 
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ti^i+i and t'2^_i— t2i-i^2i+i- To obtain the ground state, we optimize the functional 



For brevity, we set t 

E{{9i,(l)i})=Et{{0i,4)i}) +JS^J2ii j) [cos^icos^j -|-sin6'isin6'jcos(0i — 0j)] by the hybrid use of the conjugate gra- 
dient method and the simulated annealing by a monte-carlo method. Here Et{{9i,4)i}) is obtained by the numerical 
diagonahzation of (||). We should stress that the infinite system with incommensurate spin structure is embedded 
into the 2q x 2q matrix in the momentum space without finite size corrections. The numbers of the parameters to be 
optimized are reduced to Aq. 

Our recursive procedure of arrival to the ground state are the foUowings: (i) Set q — 2. (ii) Under the spin-induced 
Peierls conjecture, the ground states of x = p/q and rp/{rq) are the same, where p, q, and r are integers. We optimize 
the ground state energy for the model with g-sites in a unit cell whose lowest p sub-bands are filled, (iii) We optimize 
the ground state energy for the model with rg-sites in a unit cell whose lowest rp sub-bands are filled so that the 
difference of the energy from (ii) below a desirable precision. (We set here three or four digits.) (iv) Using the same 
procedure of (iii), we optimize the energy for p/(rg)-fillings for each q, where p and rq are coprimers. Then, we expect 
that the accuracy is the same that of (iii). (v) Setting q' to be a larger integer than q, return to (ii). The procedure 
(i-v) also provides a cross check of the spin-induced Peierls conjecture. Because this suggests the absence of a more 
larger modulated structure than g-periodic one for p/q-SWing. 

We show the numerical results of the energy spectrum in Fig. n^ where the fillings, x = p/{2q), are exhausted 
all the combinations of p = 1, 2, • • •, 17 and q = 1, 2, ■ ■ ■, 18. A Peierls gap, which arised from the spin-induced 
Peierls mechanism, opens at the Fermi surface. The structure of the energy spectrum is similar to that obtained by 
Hofstadter |18]. It is observed that the band structures are different between two regions, x < 0.25 and 0.25 < x < 0.5. 
In the former region, most of the sub-bands have enough dispersion and the energy gaps between the sub-bands are 
very small. Most of the energy gaps cannot be displayed in this magnification. In the latter, the sub-bands are less 
dispersive accompanied with larger energy gaps than those of the former. 

To see these properties in detail, the electronic properties at the irrational filling are studied as systematic limits 
of rational ones. To be specific, we focus on the cases x^^''=(/7/4 as a candidate for a;<0.25 and x^'^'>—ip/2 as that 
for 0.25<a;<0.5, where ip={^/5—l)/2. To do this, we use the Fibonacci numbers F, defined by Fn+i—Fn+Fn-i and 
Fo=Fi=l. Then the series of rational number Xn —Fn-i/ (kFn) converges to x'^^^ for A;=4 and x'-^^ for k=2 as n 
becomes large. We denote the number of the sub-bands by i3„=fcF„ {~q)- In order to specify the branching of the 
energy bands, we use sequence, {ci, C2, • • •} |16|. 

The energy spectra of the series for x^^' is shown in Fig. 0(a). The energy gaps between the sub-bands are very 
small. Except for the edge of the spectra, the width of the sub-bands become uniform as the n increases. This suggests 
that rational x=p/q gives bands with a width varying like q~^ for large q. The scaling analysis for the width, Wn, 
is shown in the upper part of Fig. ^(a). The states displayed in this figure are (1) the first sub-band just below the 



fermi surface, (or edge of the band defined by {0, 1, 1, 1, ■ • ■}) and (2) the center of the cluster shown in Fig. ||(a) (or 
center of the band, {0,0, 0,0, • • •}). It is observed that W„i?„~l (i.e. a=l) for (2), while W^„i3„~l/_B„ (i.e. a=l/2) 
for (1). This means that all these states are extended. At the band edge, a is 1/2 when the state is extended. This 
behavior conies from a remnant of van Hove singularities in Id bands. 

In contrast to those of a;'^', the spectrum for a;"' forms clusters, as n increases, accompanied with various size of 
the energy gaps. (See Fig. g(b).) The middle part of the spectra for n = 5, 6, 7 and 8 are magnified and shown in the 
inset. They show a striking similarity with those for n = 3, 4, 5 and 6. This similarity is observed as well in the other 
regions. (We should note that no self-similarity is found for x'^^'* . See Fig. ^(a).) The band width at the edges of the 
clusters seems to go down faster than q~^ . The scaling analysis for the width is shown the lower part of Fig. 0(a). The 
states displays in this figure are the 1st ({0, 1, 1, 1, • • •}), 2nd, 3rd, and 4th sub-bands from the fermi surface below. 



For the system with the structure shown in Fig 
points M . Wn is observed to decrease rapidly. T 



(b), we have to make a extrapolation in n in every two or three 
le width goes down faster than an algebraic law and is expected to 
be proportional to exp{—jBn) where the 7 is a constant. We have a = 0. Thus these states are localized. 

The results of the wave function for series for x^^^^ and 2:*^^^ are shown in Figs, ^(a) and (b) columns, respectively. 
They are the wave functions for the energy band just below the fermi surface. For x^^' , the localization length is likely 
proportional to the q, i.e. number of the sites in a unit cell. The wave function for x^'^' is localized. These properties 
are consistent with the results from the analysis of the band width. (However, the denominator of the filling is not 
large enough to make definitive conclusion only from Fig. H(a).) 

We perform the multi- fractal analysis for the wave function. The method was established in [|l9|,^. We use a 
more developed version J2^]. Extended, critical, and locahzed states are classified by amin = 1, 7^ IjO, and = 0, 
respectively. The estimates for the index a^m is shown in 0(b). It is observed that the amin for x^^' converges to 
a finite value. This means that the state is not localized. The amin for a;'^' is observed to take a small value. This 
means to the state has a tendency toward localization. These results are consistent with the analysis for the energy 
spectrum and the behavior of the wave function shown in Fig. H. 

Discussions: We have studied the energy spectra and wave function of the DE model at zero temperature. In 
studies of usual quasi-periodic systems, the quasi-periodicities are introduced by deterministic manners, such as the 
on-site or off-diagonal quasi-periodicities. In contrast to them, the self-similarity in the DE model is spontaneously 
obtained from the translation- invariant Hamiltonian P2[] . The modulation of the transfer integral is induced from 
the spin-induced Peierls instability. In the every step of the iteration toward the irrational filling, the modulation is 
expressed by q different amplitudes and, therefore, does not satisfy the Conway's theorem |23| . Therefore, the trace 
map would be described by a transcendant equation. 

The numerical results suggest the existence of two, at least, different insulating phases, a;<0.25 and 0.25<a;<0.5. 
(The scaling analysis of band width, multi-fractal analysis for the wave functions, and the features of the wave 
function and the energy spectra are mutually consistent.) These phases can be distinguished by dynamical responses, 
for example, the optical conductivity (j(uj). Roughly estimated, the optical conductivity is characterized by the the 
density of states near the Fermi level. The density of states for a;<0.25 are very different from those for 0.25<a;. For 
a::<0.25 the energy gaps between the energy bands are infinitesimally small and the density of states exist except for the 
Peierls gap. The optical conductivity is approximated by that of a usual band insulator (except for the contribution 
from the band edges of the subbands). On the other hand, for 0.25<a::, the band structure exhibits a self- similarlity. 
In the irrational limit, it is known that the optical conductivity is not well-defined in a usual sense p4[ | and is expected 
to maintain an anomalous scaling law. (To study the exponents [ p4| is very interesting. However, even if we use the 
method proposed in this article the system size and accuracy of the energy are too poor for such an analysis.) The 
monte-carlo method can be applicable for systems in d>l. In d=2, at any filling x, and for sufficiently large J, it 
was conjectured |l3] that the ground state is a flux state whose energy spectrum is similar to the present one. For 
several fillings x, the conjecture is numerically confirmed p5[ , because the phase factor in the transfer integral can 
open the Peierls gap [Eq|. However, a frustration between the incommensurate structure of the ground state and the 
periodicity of the periodic lattice arises in d>2. Due to it, to obtain the fine structures of the system are extremely 
difficult. 

The approach in this paper would be one of the fairly good starting points for introducing other degrees of freedom, 
such as the orbital |p5| , lattice, electron correlation, quantum nature of the spins, anisotropics, etc. 

The authors are grateful to Chisato Kanai, Tota Nakamura, Tatsuo Suzuki, and Kazuhiro Tada for the numerical 
correspondences, and to Yoshiko Oi Nakamura for comments on the manuscript. The authors thank the Supercom- 
puter Center, Institute for Solid State Physics, University of Tokyo for the use of the facilities. 
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FIG. 1. Energy spectrum near the Fermi surface as a function of filling x at JS^ /t = 1. Spin-induced Peierls gap opens at 
the fermi surface. We showed half bottom of the full spectrum. For x < 0.25, the energy bands are divided by small gaps. The 
gaps are far smaller than the unevenness of carbon print and they look vanishing for the eyes. 

FIG. 2. Energy spectra below the fermi surface, (a) Band structures for Xn = 1/8, 2/12, 3/20, 5/32, 8/52, 13/84, and 
21/136 are shown. This series converges to x^^K Crosses indicate the location of the band edges. The energy gaps between 
sub-bands are very small. Therefore, the respective crosses which indicate the lower edge of an certain energy band and the 
upper edge of its lower band look to overlap, (b) Band structures for Xn ~ 2/6, 3/10, 5/16, 8/26, 13/42, 21/68, and 34/110 
are shown. The series converges to x'^^K The middle clusters for 8/26, 13/42, 21/68, and 34/110 are magnified and shown in 
the inset. 

FIG. 3. (a) Plots of WnB„ against _B„ for several states. The upper three are for series x'^ . (Open squares and circles 
are, respectively, the 1st and 2nd sub-band from the fermi surface below. Crosses are for the band center.) The lower four are 
for series x^ K (Bold square, circle, triangle, and rombus are, respectively, the 1st, 2nd, 3rd, and 4th sub-bands from the fermi 
surface below.) The index of the sequence of the energy bands are {0, 1, 1, 1, 1, • • ■} for the squares and {0, 0, 0, 0, 0, • • •} for the 
crosses, (b) Estimates of the scaling index a-min for a;'^' (open square) and x'^' (bold square) as a function of 1/n, where n is 
the index of _B„. 

FIG. 4. The wave functions [^/^i [^ of the energy band just below the fermi energy, (i.e. {0, 1, 1, 1, 1, • ■ ■} series) , as a function 
of the site index i. The left and right columns show those for x*^' and x^'^\ 
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